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Abstract. We review the notion of symmetry breaking and restoration writhin 
the frame of nuclear energy density functional methods. We focus on key 
differences between wave-function- and energy-functional-based methods. In 
particular, we point to difficulties to formulate the restoration of symmetries 
within the energy functional framework. The problems tackled recently in 
connection with particle-number restoration serve as a baseline to the present 
discussion. Reaching out to angular-momentum restoration, we identify an exact 
mathematical property of the energy density E^^^ (R) that could be used to 
constrain energy density functional kernels. Consequently, we suggest possible 
routes towards a better formulation of symmetry restorations within energy 
density functional methods. 



1. Introduction 

1.1. Spontaneous symmetry breaking 

Symmetries are essential features of classical and quantal systems. For the latter in 
particular, symmetries characterize the energetics of the system and provide transition 
matrix elements of operators with specific selection rules. In nuclear systems for 
example, electromagnetic and electro-weak decays display patterns associated with 
such selection rules. 

On the other hand, certain emergent phenomena relate to the spontaneous 
breaking of those symmetries [T]. In the thermodynamic limit, i.e. when the number 
of particles N and the volume V of the system go to infinity such that N/V remains 
constant, a state with lower symmetry than the Hamiltonian can be rigorously used as 
an effective ground-state of the system. Such a state is a linear superposition of nearly- 
degenerate eigenstates, i.e. it is a wave-packet. In finite systems however, quantum 
fluctuations make such a wave-packet to relax into the symmetry-conserving ground- 
state and cannot be ignored; i.e. the concept of spontaneous symmetry breaking is only 
an intermediate description of the system that arises within certain approximations 
and symmetries must eventually be restored. Still, it makes physical sense to go 
through such an intermediate description as pseudo spontaneously-broken symmetries 
(i) relate to specific features of the inter-particle interactions, (ii) characterize internal 
correlations and (ii) leave clear fingerprints in the observed excitation spectrum of the 
system. 
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Figure 1. (Color online) Schematic representation of the different strategies 
followed by many-body methods regarding the treatment of symmetries, e.g. in 
Density Functional Theory and nuclear Energy Density Functional approaches. 



In atomic nuclei, several symmetries, if allowed, tend to break spontaneously in 
approximate descriptions based on the mean-field concept. The most important ones 
relate to the invariance of the nuclear Hamiltonian H under spatial translations and 
rotations as well as to the gauge invariance associated with particle-number symmetry. 
As described in Tab. [1] the spontaneous breaking of these three symmetries relates 
to the short-range and dominant quadrupole-quadrupole terms of the nuclcon-nucleon 
interaction as well as to its strong attraction in the L — partial-wave of relative 
motion. In particular, the strong attraction in the L — partial-wave in particular 
generates a S-wave di-neutron (di-proton) virtual state at almost zero scattering 
energy that is the precursor of neutron (proton) Cooper pairs and superfluidity in 
the nuclear medium. Even though such symmetries must be eventually enforced, 
their underlying breaking impacts the low-lying spectroscopy of finite nuclei through 
the presence of surface vibrational excitations, rotational bands and a gap in the 
individual excitations of even-even nuclei, respectively [2]. Parity and time-reversal 
are other good symmetries of H that can be spontaneously broken, while isospin 
symmetry is only approximate in the first place. 

1.2. Wave-function-based methods 

As schematically shown in Fig. [U quantum many-body methods separate into 
two categories as for the way symmetries are dealt with, i.e. (i) methods enforcing 
symmetries throughout and (ii) those that explicitly single out the intermediate 
breaking of symmetries. Although hybrid approaches that allow the breaking of some 
symmetries while enforcing the others can be set up, the present work focuses on 
those that strongly rely on the concept of symmetry breaking, i.e. methods whose 
philosophy, apart for computational constraints, is to allow all symmetries to break 



Table 1. Links between the spontaneous breaking of translational, rotational 
and particle-number symmetries and features of the nuclear force, correlations in 
the internal motion of nucleons and patterns in the excitation spectrum. 
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Figure 2. (Color online) Schematic view of the energy as a function of the phase 
and magnitude of the order parameter q of a spontaneously broken symmetry. 



spontaneously a priori. The breaking of each symmetry is monitored by the magnitude 
and the phase of an order parameter q, such that the (approximate) energy is 
independent of the phase as schematicahy shown in Fig.[2j This corresponds to the fact 
that a spontaneous symmetry breaking is accompanied by the presence of a zero-energy 
Goldstone mode. Of course, that a certain symmetry does break spontaneously usually 
depends on the number of elementary constituents of the system under consideration. 
For example, while translational symmetry (strongly) breaks in all nuclei, particle- 
number symmetry tend to (weakly) break in all but doubly-magic nuclei whereas 
rotational symmetry remains unbroken if either the neutron number or the proton 
number is "magic"|j. Figure [3] displays the correlation energy incorporated in ^"^^Pu 
and ^■^'^Sn ground-states energy through the spontaneous breaking of rotational and 
particle-number symmetries, respectively. Such symmetry breakings may account for 
up to 20 MeV correlation energy out of about 2 GeV binding energy, i.e. for about 
2%, which is much larger than the targeted accuracy on nuclear masses. Incorporating 
such correlation energies through symmetry-conserving approaches, e.g. configuration 
interaction methods, would necessitate tremendous computational efforts in such 
heavy open-shell nuclei. 

As already stated, methods authorizing the breaking of symmetries at a certain 
level of approximation must eventually restore them in a second stage. In wave- 
functions based methods, the symmetry breaking step, e.g. the symmetry unrestricted 
Hartree-Fock-Bogoliubov approximation, relies on minimizing the average value of the 
Hamiltonian for a trial wave-function that does not carry good quantum numbers, i.e. 
which mixes irreducible representations of the symmetry group of interest. Restoring 
symmetries amounts to using an enriched trial wave-function that does carry good 
quantum numbers. In terms of the schematic " mexican-hat" picture of Fig. [21 this 

I The fact that the neutron or proton number is magic is not known a priori but is based on a 
posteriori observations and experimental facts. In particular, the fact that traditional magic numbers, 
i.e. N,Z = 2, 8, 20, 28, 50, 82, 126, remain as one goes to very isospin-asymmetric nuclei is the subject 
of intense on-going experimental and theoretical investigations [3]. 
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Figure 3. (Color online) Energy gain from spontaneous symmetry breaking and 
symmetry restoration as a function of the magnitude of the order parameter q. 
Left: breaking and restoration of rotational symmetry in the ground state of 
240 ^ function of the axial quadrupole moment of the single- nucleon density 

distribution (adapted from Ref. 4). Right: breaking and restoration of neutron- 
number symmetry in the ground state of ^^''Sn as a function of the norm of the 
anomalous pair density (adapted from Ref. [5]). 



corresponds to incorporating zero-energy fluctuations associated with the phase of the 
order parainetei[|. One typical approach is to project out from the symmetry-breaking 
trial state the component that belongs to the intended irreducible representation [5]. 
Figure [3] shows that doing so for rotational and particle-number symmetries adds a 
few MeV correlation energy to the ground-state binding energy of heavy nuclei. This 
is still significant compared to the few hundreds keV targeted accuracy on nuclear 
masses. As shown in Fig. [T] a variant consists of performing the projection only 
approximately such that the calculation boils down to the minimization of a corrected 
energy expressed in terms the original symmetry-breaking state. Typical examples 
are Lipkin [HIT] of Kamlah approximate projection methods [Hl^. While it is likely 
that the strongly broken translational symmetry can be safely treated through such 
approximate projection method^, it is still unclear whether the same is true for 
weakly broken symmetries such as particle number symmetry or rotational symmetry 
in transitional nuclei. 

1.3. The nuclear energy density functional method 

Wave-function-based projection methods and their variants are well formulated 
quantum mechanically [2]. The goal of the present paper is to discuss their Energy 
Density Functional (EDF) counterparts :IT which have been empirically adapted from 
the former to deal quantitatively with properties of nuclei. The nuclear EDF method 
does rely heavily on the concept of spontaneous symmetry breaking and (approximate) 
restoration. In that sense, it is intrinsically a two-step approach. The so-called single- 
reference EDF (SR-EDF) method has been originally adapted from the symmetry- 
unrestricted Hartree-Fock-Bogoliubov method by using a density- dependent effective 
Hamilton "operator" [T^]. Later on, the approximate energy was formulated directly as 

§ Although it is not the focus of the present work, the restoration of symmetries must be 
complemented with the inclusion of collective quantum correlations associated with the fluctuations of 
the magnitude of the order parameter. This can be done through the so-called Generator Coordinate 
Method (or its energy density functional counterpart) that is formally identical to projection methods 
discussed here. 

II Such a statement is to be taken with a grain of salt for rather light nuclei |10| . 
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a possibly richer functional of one-body density matrices computed from a symmetry- 
breaking product-state of reference. The second step, carried out through the multi- 
reference (MR) extension of the SR-EDF has been adapted from the projected Hartree- 
Fock-Bogoliubov method. Such a step necessitates a prescription to extend the 
SR energy functional^ associated to a single auxiliary state of reference, i.e. a 
diagonal energy kernel, to the non-diagonal energy kernel associated with a pair 
of reference states (see Sec. I4.2p . Constraints based on physical requirements have 
been worked out that limit the number of possible prescriptions to do so [13]. Still, 
pathologies [111 [TSl US] of MR-EDF calculations have been recently identified and 
corresponding cures [TT] [181 [IS] have been proposed. Besides the actual successes 
of nuclear EDF calculations [H], the work of Refs. H?] [Hj [HI [20] demonstrates that 
nuclear SR- and MR-EDF methods must be further constrained to become satisfactory 
many-body approaches to finite Fermi systems. The first goal of the present work will 
be to reformulate, focusing on group-theory considerations, concerns about MR-EDF 
calculations that have been dealt with in Refs. [I?! [HI [IHl [20] ■ Our second objective 
will be to provide a new mathematical property that could be used in the case of 
angular-momentum restoration to constrain the form of basic EDF kernels at play. 

Given the efforts needed to better formulate EDF methods, one may question 
the necessity to stick to such approaches rather than to use a wave-function-based 
approach that strictly computes the energy from a (state-/density-independent) 
Hamiltonian, e.g. through many-body perturbation theory. Although the breaking 
(up to a few tens of MeV) and the restoration (up to a few MeV) of symmetries bring in 
both types of methods essential correlations that vary rapidly with nucleon numbers, 
incorporating the bulk of correlationf0 (hundreds of MeV) requires involved wave- 
function-based calculations (35] that are still impractical for heavy open-shell nuclei. 
The power of the EDF approach is to parameterize bulk correlations under the form 
of a functional of the one-body density (matrices) such that systematic calculations 
of heavy nuclei are tractable. The success of the overall approach, based on the 
resummation of bulk correlations into the EDF kernel and the further breaking and 
restoration of symmetries, relies on the empirical decoupling of the different categories 
of correlations at play, i.e. on the different scales that characterize them (see Tab.[2|), 
and on the fact that quickly varying correlations with the filling of nuclear shells are 
explicitly accounted for through symmetry breaking and restoratiorQ. 

Table 2. Schematic classification of correlation energies as they naturally appear 
in nuclear EDF methods. The quantity Aval denotes the number of valence 
nucleons while Gjog characterizes the degeneracy of the valence major shell. 



Correlation energy 


Treatment 


Scale 


Vary with 


Bulk 

Static collective 
Dynamical collective 


Summed into EDF kernel 
Finite order parameter q 
Fluctuations of q 


~ 8 ^ MeV 

< 25 MeV 

< 5 MeV 


A 



5 The density-dependence of the effective Hamilton operator in more standard formulations. 
+ We take as a loose definition of bulk correlations the correlation energy computed beyond a genuine 
Hartree-Fock approximation in terms of the vacuum (low- momentum I21l i nuclear Hamiltonian for 
nuclei that do not break any symmetry besides translational invariance, i.e. doubly-magic nuclei. 
* Here we have in mind to add the fluctuations of the magnitude of the order parameter. 
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1.4. Density functional theory 

As the aim of the present paper is to raise questions about the treatment of symmetries 
within the nuclear EDF method, let us make a few relevant statements about Density 
Functional Theory (DFT) [23l [2ll [25] that provides a formal framework to obtain the 
ground-state energy and one-body density of electronic many-body systems. It has 
become customary in nuclear physics to assimilate the SR-EDF method, eventually 
including corrections a la Lipkin or Kamlah, with DFT, i.e. to state that the 
Hohenberg-Kohn theorem underlays nuclear SR-EDF calculations pSl [771 1^ [5D] . 
This is a misconception as distinct strategies actually support both methods. Whereas 
the SR-EDF method minimizes the energy with respect to a symmetry-breaking trial 
density (ies), DFT relies on an energy functional whose minimum must be reached for a 
one-body density that possesses all symmetries of the actual ground-state density, i.e. 
that displays fingerprints of the symmetry quantum numbers of the underlying exact 
ground-state |31j . As a matter of fact, generating a symmetry-breaking solution is 
problematic in DFT, as it lies outside the frame of the Hohenberg-Kohn theorem, and 
is usually referred to as the symmetry dilemma. To by-pass the symmetry dilemma and 
grasp kinematical correlations associated with good symmetries, several reformulations 
of DFT have been proposed over the years, e.g. see Refs. [32l |33], some of which are 
actually close in spirit to the nuclear MR-EDF method |32j . 

Recent efforts within the nuclear community have been devoted to formulating 
a Hohenberg-Kohn-like theorem in terms of the internal density, i.e. the matter 
distribution relative to the center of mass of the self-bound system [SU [35] . Together 
with an appropriate Kohn-Sham scheme [35], it allows one to reinterpret the SR- 
EDF method as a functional of the internal density rather than as a functional 
of a translational-symmetry-breaking density. This constitutes an interesting route 
whose ultimate consequence would be to remove entirely the notion of breaking and 
restoration of symmetries from the EDF approach and make the SR formulation a 
complete many-body method, at least in principle. To reach such a point though, the 
work of Refs. [Ml |3S] must be extended, at least, to rotational and particle-number 
symmetries, knowing that translational symmetry was somewhat the easy case to deal 
with given the explicit decoupling of internal and center of mass motions. 

1.5. Outline 

The paper is organized as follows. Section [2| provides the minimal set of group 
theory considerations that are needed to go through the remaining of the paper. 
Section [3| discusses symmetry breaking and restoration within a wave-function based 
method while Sec. [4| provides the analog within the EDF context. In Sec. [4l 
we focus on difficulties encountered to formulate symmetry restorations within the 
energy functional approach. Finally, Sec. [5] briefly discusses a mathematical property 
associated with angular-momentum conservation that could used to formulate new 
constraints on off-diagonal energy functional kernels at play in MR-EDF calculations. 

2. Symmetry group 

Let us consider an arbitrary continuous compact group Q = {TZ{g)} parameterized by 
r real parameters g = {gf, i = 1, . . . , r} and whose transformations leave H invariant. 
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We denote by vg the volume of the group 

vg = dm{g) , (1) 



Jg 

where is the invariant measure on Q. Having in mind to deal more specifically 

with particle number and rotational symmetries, we further consider to be a Lie 
group, although this is not mandatory. We thus introduce the set of infinitesimal 
generators C ~ {Cf, i = 1, . . . , r} that make up the Lie algebra and in terms of which 
any transformation of the group can be expressed. The Casimir of the group built from 
the infinitesimal generators and a non-degenerate invariant bilinear form is denoted 
by A. We also denote by R{g) (C) a unitary representation of Tl{g) (C) on the Fock 
space of quantum mechanics and by S^ij{g) = R{g) the matrix elements of 

the unitary irreducible representation labeled by A. Noticing that S^ij{0) = 6ab for all 
A, the action of two successive transformations and the unitarity of the representation 
can be both read off the following identity 

E (ff) ^c.(ff) = E ^-(-ff) S^,{g) = St,{g^g') , (2) 

c c 

where —g and g—g' symbolically denote the parameters of transformations R~^{g) and 
R~^{g')R{g), respectively. States |8'*"^) span the irreducible representation A whose 
degree is dx- They are eigenstates of the Casimir A and of a chosen generator Ci 

Aie^'^) =?(A)|e^'^) , (3) 

C,\Q^^) ^ g{a)\Q^^) , (4) 

where eigenvalues 1{X) and g{a) are functions of labels A and a, respectively, with a 
running over d\ values. A so-called irreducible tensor operator and a state |8'^") 
transform according to 

R{g) R{gr^ = E St{g) , (5) 

b 

i?(3)|e'"> =Y.\Q^')Sl{g) . (6) 

b 

The discussion below is conducted for the energy, i.e. for a scalar operator H belonging 
to the trivial irreducible representation A = characterized by S^g^{g) = Sab- However, 
such a discussion can be extended to any irreducible tensor operator [36j . 

For nuclear structure, two groups are of particular importance as discussed in 
the introduction, i.e. SO{3) for rotations in the three-dimensional space and U{1) for 
rotations in the gauge space associated with particle number. The group of spatial 
translations is also essential but corresponds to a symmetry that is strongly broken 
in all nuclei and that does not need to be exactly restored in practice. Consequently, 
Tab.|3]gather useful elements that characterize U{1) and 5*0(3) such that the formulae 
given in the paper for a generic compact Lie group can be easily adapted to either of 
them. 



3. Wave- function- based method 

The present section describes what we denote as a wave-function-based method where 
energy kernels are explicitly and strictly computed as matrix elements of a Hamilton 
operator that does not depend on the wave-function it is used with, e.g. H is not an 
effective operator depending on the density of the system. 
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Table 3. Characteristics of SO{3) and (7(1) relevant to the present study. The 
gauge angle of U{1) is ip £ [0, 27r] whereas Euler angles parameterizing 50(3) 
are O = (o,/3,7) G [0,47r] X [0, tt] x [0,27r]. The one-dimensional irreducible 
representations of U{1) are labeled by m G Z whereas the (2J-|- l)-dimensional 
ones of SO{3) are labeled by 2J £ N and are given by the so-called Wigner 
functions ©^(^j^^, (O), where (2M,2M') G with -2 J < 2M, 2M' < +2J. For 
[/{!), one has 1{N) = N^, whereas for SO{3), with the choice d = Jz, one has 
1{J) = h^JjJ + 1) and g(M) = M/l. 
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3.1. Symmetry breaking 

In the present discussion, states jO"^") defined above can be eigenstates of H , which 
corresponds to handling exact solutions of the many-body problem. In general though, 
we will only consider states that are approximations to those exact eigenstates of H . 
According to our introductory discussion, a normalized, symmetry-breaking state \Q) 
takes the form 

|e) = ^c.ie^'^) , (7) 

where ^^^^ \c\a\^ — 1- Using either Eq. [5]or Eqs. |6]and[7l one can easily prove that 
the average energy 

is a scalar under all transformations of Q. However, such an energy cannot be labeled 
by good quantum numbers (A, a), which is the fingerprint of the symmetry-breaking 
character of the many-body state |0). 

3.2. Symmetry restoration 

One extends the definition of the symmetry breaking energy E (Eq. |5]) to the so-called 
energy kernel E[g' , g] through 

{e\R-Hg')HR{g)\Q) 

Hg.g\- (e|i?-i(5')i?(5)|e) ' 

where the norm overlap kernel is N[g',g] = {Q\R^^{g')R{g)\<d) . Expanding |8) 
according to Eq. [7] and using Eq. [2l one obtains 

E[g',g]N[g',g] = Y^^la^xb S^,{g-g') , (10) 

Xab 

N[g\g] = , (H) 

Xab 

where the symmetry-restored energies {Q^"'\H\Q'^ " ) = E^ 6\x> Saa' provide the eigen- 
spectrum of H if one is working with its exact eigenstates. Expressions |9] and [TT] 
correspond to the double expansion over the volume of the group 

^[5'' 5] ^ E E ^abl'b' S^,{g') S^,,, (5) , (12) 

AA' aba'b' 
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applied to functions F[g',g] — f[g — g'] that in fact only depend on the difference of 
the two arguments. In such a case, the double expansion reduces to a single expansion 
whose coefficients /^^ are e^^ = c^^ cxf, and = cxb for the functions of 
interest e[g—g'] n[g—g'] and n[g—g' .respectively. Given that coefficients /^j, transform 
specifically under any member of Qfl, the ratio of two such objects, e.g. E^, transforms 
as a scalar. 

Starting from E[g',g] and N[g',g], and given the orthogonality relationship 

dmig) Sab (g) Sa'b' (.9) = ' (1^) 

g "A 

one can perform the integration 

^) JJ/'^(9')dmig)Si{g')Si*ig)E[g',g]N[9',g]^e'a, , (14) 

to extract the energy E^ associated with states jO'^") spanning the irreducible 
representation Such a symmetry restoration stage is denoted as a multi-reference 
method in the sense that, while the energy computed through Eq. |S] involves a 
single reference state |0), the extraction of E^ involves the set of references states 
|9(g)) = R{g)\Q) obtained from \Q) through all transformations of Q. 

3.3. Transfer operator 

The above presentation of the symmetry restoration will be particularly suited to 
making the connection with the MR-EDF approach discussed in Sec. ID Still, it is 
worth noting that Eq.[T4]is usually obtained within the wave- function-based approach 
by first introducing the transfer operator [2] 

P^,^^ [ dm{g)Sl:{g)R{g) , (15) 
vg Jg 

to extract the state with good symmetries 

|0'°> = — ^a'j0) , (16) 
C\b 

from which Eq. [T3]can be recovered by taking the average value of H. 

3.4. Particular case of interest 

Let us consider the particular case where the symmetry-breaking state \Q) is taken 
as a product state |$) of the Bogoliubov type. In this case, the method at 
play corresponds to symmetry-unrestricted and symmetry-projected Hartree-Fock- 
Bogoliubov approximations [5]. In such a situation, one can use the standard Wick 
theorem [Sj to express the symmetry-breaking energy E (Eq.[8]) as a specific functional 
E[p, K, K*] of the one-body density matrices computed from the symmetry breaking 
state and given, in an arbitrary single-particle basis, by 

_ (<i>|a]a,|<f>) _ma,a,\<i>) , _ m4a]m 

'^a = — — ■ U'J 



- ($1$) ' " ($1$) ' ~ ($1$) 

(t The corresponding law is easily obtained from the transformation of S^f^{g). 

ffThe fact that E[g',g] and N[g',g] only depend on the difference g — g' can be exploited to extract 
e^jj through a single integral rather than through a double integral as in Eq. 1141 The reason why we 
keep explicitly two integrals here will only become clear in the last section of the paper (see Sec. 15. id . 
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The explicit form of E[p, k, k*] depends on the specific Hamiltonian _ff , e.g. on whether 
it is local or non local in space and on whether it contains a three-body force, a four- 
body force. . . 

As for the symmetry restoration stage, the energy kernel at play (Eq. [TT|) 
can be computed in this case using the generalized Wick theorem [38j such that 
E[g',g] = E[p^ ^ , ^ , *], i.e. the transition energy kernel is expressed through 
the same functional as the symmetry-breaking energy except that diagonal symmetry- 
breaking one-body density matrices are replaced by transition one-body density 
matrices involving the two transformed product states \^{g)) and \^{g')), i.e. 

g'g _ m9')Whm9)) ( $ (ffQ | | $ (ff ) ) ($ (ffQ | ^1 fl] | $ (ff)) 



3.5. Actual implementation 

In practice, states jO"^") only provide approximations to exact eigenstates. This is for 
instance the case when starting from a symmetry-breaking product state as discussed 
in Sec. 13.41 When dealing with a non abelian group, one must actually consider 
an arbitrary linear combination of states spanning a given irreducible representation 
such that mixing coefficients are determined through the minimization of the resulting 
energy. This corresponds to considering that the link between the symmetry-restored 
states of interest and the symmetry-breaking one is in fact given by 

|0Aa^^^^A6pA|Q^ , (19) 
b 

rather than by Eq. 1161 and to determining the {g^''} through the minimization of 

^ (e^"|g|9^°) ^ Ebb' g''* 9''' el. 
- {e^^\e^-) Y.bb'9'''*9^''<' ' 

with e^i defined by Eq . [T4l and n^, given by a similar equation for the kernel N[g' , g]. 
Such a minimization leads to solving a Hill- Wheeler- Griffin equation [Ml HO]. 



4. EDF-based method 



We are now interested in the nuclear EDF formalism within which energy kernels are 
formulated as functionals of the one-body (transition) densities defined by Eqs. [T7] 
and [IHl Although extensions can be envisioned, the standard EDF method is 
based on the use of simple symmetry-breaking product states of the Bogoliubov 
type. This choice stems from the historical construction of the nuclear EDF 
method by analogy with symmetry-unrestricted and symmetry-projected Hartree- 
Fock-Bogoliubov approximations fSec. 13. 4|) 



4.1. SR- EDF step 

The SR-EDF method [TT] relies on computing the analog to the symmetry-breaking 
average energy E (Eq. |51) as an a priori general functional £[p,K,K*]. As opposed 
to what was considered in Sec. 13. 4[ the symmetry-breaking energy £[p,K,K*] is not 
computed from the average value of a genuine operator H. As a result, specific 
constraints must be imposed onto the functional form of k, k*] for it to be a scalar 
under all transformations of G', i-e. under transforming |$) and the densities p, k, 
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K* constructed from it. For a quasi-local functional of the Skyrme type, we refer 
the reader to Refs. [lU im for the formulation of such constraints. For actual 
parameterizations of the nuclear EDF, we refer the reader to Ref. [TT]. 

4.2. MR-EDF step 

The MR-EDF step is built from the SR-EDF by analogy to the wave-function-based 
method formulated in terms of product states fSec. l3.4| ). As a result, the analog of the 
energy kernel E[g' (Eq. 1111) is naturally introduced over the volume of Q through 
the definition E[g',g\ = f [pf's, K^'f , k^s' *], where £[p,Hi,K*] is the SR-EDF kernel. 
From £[g' ,g\^ one extracts 

j dm{g')dm{g)Sl{g')S^C{9)£[9\9]NW,g] , (21) 

VgJ Jgjg 

by analogy to Eq.[T31 Whereas in the wave- function-based method one could explicitly 
demonstrate the identity e^^ — c\b , this is not the case in the EDF approach 
where there no possibility in general to perform the equivalent to the derivation that 
led to Eq. [9l Equation [21] corresponds simply to the application of expansion [12] to 
the function £[g',g] over the irreducible representations of the group, without any 
reference to many-body states with good quantum numbers. As a matter of fact, and 
contrarily to what is often stated jllj . e^;, is not computed from a projected state in 
the MR-EDF method, i.e. the transfer operator P^^ cannot be factorized explicitly in 
Eq. [5T] However, one can implicitly relate the MR-EDF energy £^ to the projected 
state \^^°') obtained from |$) as in Eq. [T^] With this in mind, it is natural and 
customary [111 117) to define the symmetry-restored energy £^ from e^^ through the 
analog of Eq. [20] i.e. 

f ^ ^ ^ 5''* el, I g'"* n^, , (22) 

bb' bb' 

where the {g^^} are determined through the minimization of 
4-. 3. Puzzling questions 

We have clarified in previous sections that SR- and MR-EDF methods have been 
empirically constructed by analogy to symmetry-unrestricted and symmetry-projected 
Hartree-Fock-Bogoliubov approximations. The key difference with the latter is that 
the energy kernels at play in the EDF method are not defined as matrix elements of 
a genuine operator between wave-functions. For the rest, expressions utilized in both 
approaches, in particular regarding the extraction of the symmetry-restored energy 
fEqs. [T4l and [20l versus Egs. [2T] and [22]) . look totally alike. Still, puzzling questions 
remain to be raised. 

As mentioned in Sec. 14.11 one must require at the symmetry-breaking level 
that the SR energy £[p,K,K*] is a scalar under all transformations of Q. Such a 
requirement has led to formulating a set of constraints on the functional form of 
£[p,K,K*] [41] [42] [43]. The next question one may ask is the following: are those 
constraints imposed on the energy kernel £[p, k, k*] at the SR level sufficient to making 
the MR-EDF method described in Sec. l4.2l well defined, in particular from a symmetry 
standpoint? In particular, one may wonder whether the fact that the energy kernel 
£[g',g], which is the key ingredient to the MR-EDF calculation, is not computed as 
the matrix element of a (genuine) operator makes the method ill-defined in any way? 
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As a matter of fact, a set of physical constraints to be imposed on £[g',g] have 
already been worked out [13]. The facts (i) that the MR energy should be real, (ii) 
that the SR-EDF should be recovered from the Kamlah expansion and (iii) that the 
Random Phase Approximation based on the SR-EDF £[p, k, k*] should be recovered as 
a limit of the MR-EDF calculation has helped limiting the energy kernel £[g',g] 
to depend on transition densities only, e.g. £[g\g] = £[p^ ^ , ^,k^^ *]. 

The aim of the present contribution is to elaborate further on the question raised 
above and to discuss a path that could be followed to constrain more tightly MR-EDF 
calculations. References [111 HZl [THl [12] have already provided important elements in 
the case of C/(l), i.e. for particle- number restoration (PNR). Let us recall the main 
outcome of those studies prior to formulating the problem to be addressed. 

4- 4- Lessons learnt from particle-number restoration 

Equation [11] applied to U{1) provides the Fourier decomposition e[iy9]n[(p] = 
J2n£Z '^'n e*^''' of the periodic function £[ip\ n[ip\ over [0, 2tt\. From a mathematical 
standpoint, the sum runs a priori over all irreducible representations of the group, i.e. 
over both positive and negative integers N . From a physics point of view though, the 
label N denotes the particle number of the physical system. Consequently, the sum 
should actually only run over positive integers, i.e. one should find cj^S^ = and 
£^ = ior N < 0. In the wave-function-based method, such a result is indeed obtained 
from the fact that is computed as the average value of H in |$^), the latter being 
zero [H] for iV < 0. In the EDF context, however, it was demonstrated [HI [1^ that 
Fourier components £^ are a priori different from zero for N < 0, i.e. one usually 
obtains a non-zero symmetry-restored energy for negative particle numbers! This 
problem was shown [18 to be related to unphysical mathematical properties of £ [ip] . 
Applying the regularization method proposed in Ref. |17) , the cancelation of non- 
physical Fourier components was recovered [18]. At the same time, components £^ 
for N > were modified by up to 1 MeV, which is of the same order as the root-mean- 
square error on mass residuals reached by the best available particle-number-restored 
EDF mass fits [IS]. This demonstrates the practical need of constraining further 
MR-EDF calculations in order to produce fully reliable results. 

5. Towards new constraints? 

The example discussed above is particularly enlightening given that clear-cut physical 
arguments can be used to argue that certain coefficients in the Fourier expansion 
e[tp] n[ip] should be strictly zero, although they are not if one does not pay particular 
attention to it. Recovering such physical features removes at the same time non- 
physical contaminations from other coefficients of the expansion 18 . This proves 
that the MR-EDF method, as performed so far, faces the danger to be ill-defined and 
that new constraints on the energy kernel £[g' , g] must be worked out in order to make 
the method physically sound. The regularization method proposed in Ref. |17| that 
restores the validity of PNR can only be applied if the underlying EDF £[p,K,K*] 
depends on integer powers of the density matrices [19] , which is an example of such a 
constraint. 

For an arbitrary symmetry group, the situation might not be as transparent as 
for C/(l). Indeed, it is unlikely in general that certain coefficients of the expansion of 
£[g' , g]N[g' , g] over irreducible representations of the group are zero based on physical 
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arguments. The challenge we face can be formulated in the following way: although 
expansion [T^ that underlines the MR-EDF method is sound from a group-theory point 
of view, mathematical properties deduced from a wave-function-based method must be 
worked out and imposed on £[g' , g] to make e^j, extracted from Eg . 12 1 I phvsicallv sound. 
The rest of the present contribution is dedicated to briefly introducing an example of 
such a property in the case of SO{3), i.e. for angular momentum restoration, that 
could be used to constrain 5[f2',ri]. Details of such an analysis will be reported 
elsewhere l36l . 



5.1. Mathematical property associated with angular-momentum conservation 

We omit spin and isospin for simplicity and consider the rotationally-invariant nuclear 
Hamiltonian H = T + V in which the central two-nucleon interaction 

V=^JJ dridf2 vilr'i - r2|) 4i "-k ap^ > (23) 

is local, i.e. non-antisymmetrized matrix elements are defined as (1 : ri; 2 : r2|^|l ^ 
rij; 2 : fi) = v{\fi — r2|) Si^fi -— f^) d{f2 — fi), and in which three-nucleon and higher 
many-body forces are disregarded for simplicity. None of the conclusions drawn below 
would be modified by the inclusion of many-body forces or by using a non-local two- 
nucleon interaction. Operator at [apj creates (annihilates) a nucleon at position 
f. Considering an eigenstate of and L^, as well as using center of mass 

R= {fi -|- r2)/2 and relative coordinates f=ri — r2, the potential energy reads as 



^eLAf|0LM^ 2 J J " Rf 

dRV^^'iR) , (25) 



which defines a potential energy density V {R) in terms of the two-body density 
matrix p^l,^*"*^ = (e^^^| ot at a^., a^J 9^*0/(0^*^ I ©^^0- After tedious but 



straightforward calculations, one can demonstrate that 



2L 



V^'^iR) ^ J2 ^lMl'o YL'iR) , (26) 

where the Clebsch-Gordan coefficient C'f'l^/^.'o carries the dependence on M while YJ™ 
denotes spherical harmonics. The weight u^|^,(i?) depends on the norm of R only 
and is related to a reduced matrix element of the two-body density matrix operator 
recoupled to a total angular momentum L' . The remarkable mathematical property 
identified through Eq. [55] is that the scalar potential energy V'" is obtained from 
an intermediate energy density V^'^^R) whose dependence on the orientation of R 
is tightly constrained by the angular-momentum quantum number of the underlying 
many-body state |8^*^), i.e. its expansion over spherical harmonics is limited to 
L' < 2L. Such a result is unchanged when adding the kinetic energy (density) to the 
potential energy (density) such that we restrict ourselves to the latter for simplicity. 
Of course, the energy eventually extracts the coefficient of the lowest harmonic, i.e. 
^ V^JdRv^lliR). 
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5.2. Wave-function-based symmetry-restoration method 

Since property [55] is general, it can also be obtained within the frame of the wave- 
function-based symmetry-restoration method presented in Sec. [31 Omitting again the 
kinetic energy for simplicity and using Eqs. [2]|7l lTTIl23l the potential energy kernel 
reads 

v[n',n]N[n',n] = ]- ( dRdfv{r) {e\Ri {n') p^^lR{n)\e) (27) 

= \T. cl^N,CL,^.D'^,%^{n')Dll^^{n) j dRdrV{r)p%'^^''^'^^^'^ , 

{LM} 

where {L,M} denotes a sum over the six angular-momentum quantum numbers 
appearing in the formula. Applying Eq. [T4]to the above expression (Eg. I27p provides, 
thanks to the orthogonality property [131 the result 

\2 r nL /n'\ nL* 



= + / d^i'dn £j<Mi^ RkmM. v[n', n] N[n', n] 

^\jdRdrV{r)p%^''''^^ , (28) 

so that Eqs. [25l and [26l are recovered. To obtain such a result it is mandatory to 
use the double-integral formulation of Eq. [Ml rather than the more standard single- 
integral formulation that takes advantage, from the outset, of the fact that V^[r2',r2] 
and N[Q,' , il] only depend on the difference VI— U! . We thus insist on using the double- 
integral formulation in the present discussion. 



5.3. EDF-based symmetry-restoration method 

Let us now come back to the EDF method described in Sec. [H The point is to 
underline the fact that property [26l cannot be derived a priori given that the potential 
energy part of the kernel is not explicitly related to the two-body density 

matrix in this case. Taking a quasi-local Skyrme EDF as an example, although this 
can be easily adapted to non-local EDF of the Gogny type, the energy kernel takes 
the form 

= J dR£{p^'"{R),T^'"{R),j^'^{R),...) , (29) 

where £ {p^' ^ (R) , t^' ^ (R) , j^^' {R)) is a general function of a set of one-body local 
transition densities built from the transition density matrices, e.g. 

p"'"(i?)^ ^^*(i?)^,(i?)pg.'" , (30) 

r"'"(i?) ^ J2 ■ [^^^(R)] f^'' ' (31) 

r^^iR) ^ - ^ E {^.*(^) [^^'^i^)] - [V^*(i?)] V^{R)} Pg'" ,(32) 

ij 
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such that constraints imposed at the SR level [HI 33] are fulfilled (see Sees. 14.11 
and 14. 3p . Given such an EDF, there is no reason a priori that the energy density 
£^^^{R) extracted from Eqs. [H] and [25] displays property [25] i.e. the angular 
dependence of is likely to display harmonics Yl,{R) with L' > 2L. One 

might argue that it is not an issue considering that the symmetry-restored energy £^ 
eventually relates to the harmonic Yq(R) only. However, a formalism that provides 
f^*^(_R) with a spurious angular content will certainly also provide the coefficient 
£lo{R) of the lowest harmonic with unphysical contributions. To state it differently, 
it is likely that constraining the MR-EDF kernels £[p^ ^, ^, *] to produce an 
energy density £^^\R) that fulfils the mathematical property [26l will impact at the 
same time the value of the weight £lo{R), and thus the value of £^. To some extent, 
this is similar to the situation encountered with U{1) where restoring the mathematical 
property that Fourier coefficients £^ with TV < should be strictly zero impacted the 
value of all other Fourier coefficients [H] . 

6. Conclusions 

The present contribution reviews the notion of symmetry breaking and restoration 
within the frame of nuclear energy density functional (EDF) methods. Multi-reference 
(MR) EDF calculations are nowadays routinely applied with the aim of including 
long-range correlations associated with large- amplitude collective motions that are 
difficult to incorporate in a more traditional single- reference (SR), i.e. "mean-field", 
EDF formalism [IT]. 

The framework for MR-EDF calculations was originally set-up by analogy with 
projection techniques and the Generator Coordinate Method (GCM), which are 
rigorously formulated only within a Hamiltonian/wave- function-based formalism [2]. 
In the present work, we elaborate on key differences between wave- function- and 
energy-functional-based methods. In particular, we point to difficulties encountered 
to formulate symmetry restoration within the energy functional approach. The 
analysis performed in Ref. [18] to tackle problems encountered in Refs. [14l[15][l6] for 
particle number restoration serves as a baseline. Reaching out to angular-momentum 
restoration, we identify in a wave-function-based framework a mathematical property 
of the energy density E^^^ (i?) associated with angular momentum conservation that 
could be used to constrain EDF kernels. Consequently, possible future routes to 
better formulate symmetry restorations in EDF-based methods could encompass the 
following points. 

(i) The fingerprints left on the energy density E^^^ [R) by angular momentum 
conservation in a wave-function-based method could be exploited to constrain 
the functional form of the basic EDF 

(ii) The regularization method proposed in Ref. [T7] to deal with specific spurious 
features of MR-EDF calculations should be investigated as to what impact it has 
on properties of the energy density £^^\R) in the case of angular momentum 
restoration. 

(iii) Similar mathematical properties extracted from a wave-function-based method 
could be worked out for other symmetry groups of interest and used to constrain 
EDF kernels. 

Efforts in those directions are currently being made |36j . 
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